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Abstract 

Aims. We present the results of a pilot project designed to study the distribution of dark matter haloes out to very large radii in spiral 
galaxies. As dynamical probe we use their rotation curves and the motions population of satellite galaxies. In this pilot stage, we 
observed seven late-type spiral galaxies of about the same luminosity Mr ~ —22 (and approximately the same mass). We investigate 
the kinematics of these galaxies, and the radial and angular distribution of their satellites. 

Methods. Using VIMOS at the VLT, we carried out a spectroscopic survey in seven 14' x 14' fields each around a late-type isolated 
spiral galaxy. We obtained radial velocities and spatial distributions for 77 candidate satellites. After removing the interlopers, we are 
left with 61 true satellites. In combination with the rotation curves of the primary galaxies, satellites are used to probe the gravitational 
potential of the primaries and derive the dark matter halo properties by means of standard mass modeling techniques. 
Results. We find (a) that the dark matter haloes of luminous spirals (Mr ~ —22) have virial radii of ~ 400 kpc and virial masses of 
3.5 x 10 12 M Q ; (b) that the radial velocity and angular distributions of the satellites around the primaries are isotropic; and (c) that 
the resulting mass distribution is in good agreement with that found in the optical regions of spirals and described by the universal 
rotation curve of spirals once extrapolated to large radii. The results obtained in this pilot phase of the project are already interesting 
and limited only by small number statistics. The full project involving an order of magnitude more targets, would very likely provide 
us with a definitive picture of the dark matter distribution around spirals out to their virial radii and beyond. 

Key words. Galaxies: kinematics and dynamics - galaxies: haloes - galaxies: dwarf, galaxies: structure - cosmology: dark matter 



1. Introduction 

In the past few decades, it has been revealed that dark mat- 
ter (DM) haloes extend far beyond the optical boundaries of 
galaxies. However, their actual "sizes" and, above all, their den- 
sity distribution outside the luminous regions, are still quite un- 
known. The inner regions of DM haloes are studied using several 
methods including the analysis of both rotation curves for spi- 
ral galaxies (Rubin et al. 119801 Bosma [19811 Persic, Salucci & 
Stel 1996 ), and velocity dispersions for ellipticals (Kronawitter 
et al. EOOQl Cappellari et al. 120061 Thomas et al. l2QTj7l ). Different 
tracers of the gravitational field, must be used to investigate the 
external regions, which are relatively void of luminous matter. 
These include planetary nebulae (Coccato et al. 2009) and glob- 
ular clusters in ellipticals (Spitler et al. 2009), weakly lensed 
background galaxies, and the satellites populating the outer parts 
of the DM haloes. These objects are likely to be the remnants 
of the halo formation process and provide us with an efficient 
dynamical probe, out to the "edge" of the halo and beyond, as 
they are bright, numerous, and radially extended test particles. 
By measuring the satellite-primary relative velocity for a suf- 
ficient number of objects, we can derive the mass distribution 
of the halo out to very large radii. The results are complemen- 
tary to those obtained from the analysis of weak lensing data, 
not only because the measurements involved and their related 
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biases are independent, but because the satellite motions lead to 
a measure of V(R) — (GM/,(/?)//?) 5 i.e. the halo mass inside 
a given radius, while the observed tangential shear derived the 
weak lensing measures a more complex quantity, the (average - 
local) halo surface density at a given radius £(< R) - Since 
{GM h (R)/R)°- 5 in the regions under study at most decreases as 
R~ 25 , the sensitivity of the dynamical measurements performed 
using satellite kinematics does not depend on radius, in contrast 
to what occurs for the majority of tracers of the gravitational po- 
tential including weak lensing. Finally, comparing different mass 
models derived from various kind of measurements, provides a 
valuable consistency check. 

The virial radius /?,,,> is an important characteristic scale- 
length of dark matter haloes. In the current concordance cos- 
mology (where £Icdm ~ 0.3, ~ 0.7, and h — 0.72 represent 
the dark matter, dark energy densities, and Hubble parameter, re- 
spectively), is defined as the radius that encompasses 330 times 
the background mass 4 /2>nR? vil SlcDuPc, with p c ~ 10~ 29 g/an i 
the critical density of the Universe. For historical reasons /?200, 
i.e. the radius that encompasses a mass 4/37r/? 3 , iV 200p c , is also 
often considered, where is related to Z? VI > all halo structural prop- 
erties (concentration and mass) and can be taken as the refer- 
ence size of a DM halo (e.g. Navarro et al. 120 101 ). Moreover, 
simulations show that it approximately marks the transition be- 
tween the virialized matter of a halo and the infalling mate- 
rial. In detail, for a virialized structure of mass M V! >, at red- 
shift zero, R vir = 26O(M Wr /(lO 12 M )) 1 / 3 kpc defines the typical 
scale-length relevant to this study. 
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From an observational point of view, the distribution of the 
DM outside the inner luminous regions out to /?,.,, is still poorly 
known, and data suitable for investigating these issues are scarce. 
Invaluable data is provided by weak lensing measurements (e.g. 
Mandelbaum 2008 ), which indicate that /?,,,> is some hundreds 
of kpc, although the precise value is uncertain because of the 
weakness of the signal at these distances coupled with the un- 
certainties in the DM density profile. 

In this pilot study, we show that it is possible to investigate 
the DM distribution of spiral galaxies out to R v j r by combining 
the internal kinematics of the spirals and the radial velocities of 
their satellites. 

It is well known that the spatial and kinematical distribu- 
tion of satellites around their host galaxies provide important 
information about the DM halo properties. Previous investiga- 
tions have also noted that further analysis is possible only if the 
angular distribution of satellites is homogeneous. It is extremely 
difficult to recover the gravitational potential of the galaxy halo 
from the kinematics of satellite systems if their distribution is 
not (roughly) isotropic. Zaritsky et al. d 19971 ) studied a sam- 
ple of 69 spiral galaxies and their 115 satellites, and found a 
tendency for the satellites to lie preferentially close to the mi- 
nor axis of their hosts, a phenomenon called the Holmberg ef- 
fect (Holmberg [19691 . Brainerd (120051 hereafter B05) studied 
a sample of 2000 primaries and their 3300 satellites selected 
from the Sloan Digital Sky Survey Data Release 3 (SDSS DR3; 
York et al. 2000), and found instead a tendency for satellites to 
lie near the major axis, a result confirmed by Faltenbacher et 
al. d2007b using SDSS DR4 data. Agustsson & Brainerd d20"l0l ) 
studied a large sample of SDSS DR7 galaxies (4487 hosts and 
7399 satellites) and for blue hosts found that at projected radii 
r p < 150 kpc satellites are distributed close to the major axes, 
while at r p > 300 kpc satellites are found close to the minor axes. 
Azzaro et al. (2006) found an isotropic distribution of satellites 
based on a sample of 193 satellites and 144 isolated spiral galax- 
ies, though they issued a caution about the insufficient quality 
of their statistics. Using SDSS DR5 data, Bailin et al. ([20081 
confirmed this result with 273 late-type primaries hosting 321 
satellites. 

We believe that the contradictory results found so far have 
the following causes. First, most of these studies rely on data for 
a very small number of satellites per host. The average primary 
galaxy often possesses data for only two satellites, which is a 
very small fraction of the true total number, as we can see just 
by looking at the Local Group (Hartwick 2000; Kroupa et al. 
[20051 Willman et al. 120031 Koch et al. 120061 . Second, to com- 
bine the kinematics of systems composed of primaries of differ- 
ent masses and halo sizes is very complex and the subsequent 
analysis is far from trivial and can become easily biased. 

This uncertain situation has suggested to us a new observa- 
tional approach to the problem. First, we select a relatively small 
number of primary galaxies (compared to previous studies), and 
all in a quite narrow range of luminosity —21.6 < Mr < —22.7. 
For each of them, we then detect a larger number of satellites 
by using powerful instruments (VIMOS at the VLT). The char- 
acteristic halo mass will eventually be sampled by 30 primaries 
in that luminosity range, and their gravitational potential will be 
reconstructed by means of a) the disk kinematics and b) the kine- 
matics of their about 300 possibly detected satellites. 

In this pilot paper, which is also a feasibility study, we carry 
out about 25% of the project. To complete the project and inves- 
tigate the mass distribution around spirals of different masses, 
we plan to replicate the investigation around objects of signifi- 
cantly lower luminosity (e.g. with —20 < Mr < —20.5). 



We measure the rotation curves (RCs) of seven primaries to 
help us reconstruct the mass distribution in the inner regions of 
galaxies and eventually ensure that the systems have approxi- 
mately the same total mass (2 — 5 x 10 12 M©). We then measure 
the primary-satellite relative velocities for the candidate satel- 
lites and, after removing interlopers, we build for the latter the 
velocity dispersion profile of the co-added system. The method 
has a number of advantages: a) Since the primaries have almost 
the same mass, we can co-add the kinematical measurements 
from different objects as they would belong to a single one. This 
leads to a very precise measure of the rotation curve in the opti- 
cal region, a good hint of the satellite number profile, and a fair 
measure of their velocity dispersions, b) Since we detect objects 
down to a very low (Mr = — 16) limiting magnitude, our dynam- 
ical analysis is not affected by possible luminosity-dependent bi- 
ases in the positions and in the motions of satellites, c) We can 
easily reach the virial radius and beyond. 

In this paper, we report the first results of the pilot phase 
of our project, obtained with seven hours of VIMOS at the 
VLT, followed later by a program totaling 53 hours at the VLT 
(VIMOS) and NTT (EMMI). In this phase, we obtain valuable 
information about the properties at the outer radii of dark and lu- 
minous components around spirals, but we are also particularly 
interested in demonstrating the feasibility of the project. We an- 
ticipate that, once completed, our project will be able to unravel 
most of the distribution of DM within the haloes surrounding 
spiral galaxies out to their "edges". 

The layout of the paper is the following. In Section 2, we de- 
scribe observations and data reduction of primary galaxies and 
their satellite candidates. In Section 3, we perform the data anal- 
ysis: in Section 3.1 of the structural properties of the primaries, 
using their rotation curves; in Section 3.2 identifying interlop- 
ers in the satellite candidates; in Section 3.3 studying the radial 
distribution of satellites; in Section 3.4 studying the angular dis- 
tribution of the satellites; in Section 3.5 estimating the velocity 
dispersion profile; and in Section 3.6 performing mass modeling 
of the primary galaxies. We discuss the results in Section 4. 

2. Observations and data reduction 

The selection criteria of the primary galaxies were the following: 
(1) all galaxies are late-type spirals; (2) they have an inclination 
between 40° and 70°; (3) they are isolated: all objects within a 
radius of 700 kpc and 1000 kins" 1 in relative radial velocity 
have a luminosity at least 2.5 magnitudes below that of the pri- 
mary galaxy; (4) they have a redshift 0.029 < z < 0.085; and (5) 
they have at least four known satellites (from the SDSS). To se- 
lect the isolated primaries, we used the catalog of host galaxies 
and satellites from B05. 

The observations were conducted with two instruments: the 
kinematics of primary galaxies were measured using data ac- 
quired with EMMI at the NTT, while the hunt for satellites 
and the measurement of their radial velocity was done with the 
VIMOS at the VLT. Redshifts and luminosities of the primary 
galaxies, and the observing log, are presented in Table Q] 

2.1. Kinematics of primary galaxies 

The major-axis ionized gas rotation curves of the primary galax- 
ies were obtained by means of long-slit spectroscopy at the ESO- 
NTT 3.5m telescope in La Silla (077.B-0767), using the ESO 
Multi-Mode Instrument (EMMI) in the REd Medium Dispersion 
(REMD) configuration. Grating No. 6 with 1200 grooves mm -1 
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Primary galaxy z mag a Exp.Time EMMI b Exp.Time VIMOS b 



J003 828+000720 


0.042 


14.60 


2 x 2200 


2 


X 


2180 


J134215+015126 


0.029 


13.60 


2 x 2000 


5 


X 


2180 


J14521 1+044053 


0.046 


14.62 


3 x 2000 


6 


X 


2180 


J 15262 1+035002 


0.083 


14.94 


2 x 2000 


5 


X 


2180 


J153221-002549 


0.085 


15.18 


2 x 2200 


8 


X 


2180 


J154040.5-000933.5 


0.076 


15.70 


4x2700+1800 


6 


X 


1000 


J 154904-004023 


0.077 


15.73 


2 x 2200 


5 


X 


2180 


J221957-073958 


0.038 


13.57 


1800 


5 


X 


2180 



a r-band magnitudes from SDSS 
k exposure time is in seconds 



* * , V** 
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Figure 1. Rotation curves of the seven primary galaxies depro- 
jected for the disk inclination and folded around their centers. 



blazed at 6890 A was used in the first order in combination 
with the two red-arm mosaicked MIT/LL and a l"x330"slit. 
Data reduction was done with MIDASEJand IRAlS standard rou- 
tines. Wavelength calibration was performed using observations 
of a helium-argon calibration lamp. All frames were bias sub- 
tracted, flat-field corrected, cosmic-ray cleaned, and wavelength 
calibrated. Quartz lamp and twilight sky flat-fields were used 
to remove pixel-to-pixel variations and large-scale illumination 
patterns. After wavelength calibration, repeated exposures of the 
same target were aligned and summed. We used the night-sky 
emission lines to assess the quality of the wavelength calibration 
and spectral resolution. We use these data not only to derive the 
rotation curve of the galaxies but also to define their recession 
velocity. This velocity is later compared with the one measured 
for the satellites using VIMOS. The instrumental resolution is 
1 .3 A (FWHM), corresponding to a resolution of a = 27 km s~ 1 
at Ha. The ionized gas kinematics was measured by simulta- 
neous Gaussian fits to Ha, [Nil], and [S II] lines by means of 
a self-written IDL routine and then folded around their centers 
and deprojected for the disk inclination as done by Pizzella et al. 
(2004). The derived velocity curves are plotted in Fig.Q] 



MIDAS is developed and maintained by ESO. 
2 IRAF is distributed by the National Optical Astronomy 
Observatories which are operated by the Association of Universities for 
Research in Astronomy under cooperative agreement with the National 
Science Foundation. 



2.2. Satellites radial velocity 

Observations were carried out in service mode with the Visual 
MultiObject Spectrograph (VIMOS) operated at Cerro Paranal 
(075.B-0794; 077.B-0767). To identify satellites orbiting se- 
lected primary galaxies and measure their radial velocity, the 
instrument was used in MOS configuration. After pre-imaging, 
we constructed the masks for the MOS observations selecting 
the brightest non-stellar objects in the field. 

We adopted grism "HR orange" on the four arms of the in- 
strument. The grism is used in combination with the order sort- 
ing filter GG435. The grism is characterized by a dispersion of 
0.6 A px~', which corresponds to a spectral resolution of 2150 
for a 1" slit. We could accommodate, on average, about 50 slits 
in each quadrant for a total of about 200 slits for each primary 
galaxy. We organized the observations by repeating the same 
pointing with the same mask several times. Each exposure is 
2180 seconds long. Some exposures that did not fulfill our ob- 
servational requirements were repeated. For this reason, we ob- 
tained more exposures than requested. Extra exposures, executed 
under weather conditions not complying with our constraints, 
were only used if the data quality was good. The exposures re- 
ported in Table[T]are the ones used in our final data reduction. 

We started our data reduction from frames that were bias 
subtracted, flat-field corrected, and wavelength calibrated by the 
ESO data reduction pipeline. Further data reduction was done 
with MIDAS. The first step was to assess the quality of each inte- 
gration. Once the selection was done and the low quality spectra 
discarded, we computed the average spectrum for all frames. The 
second step was to approximately derive, in the sky-subtracted 
spectra, the redshift of all targets. To this aim, we used emission 
and absorption lines listed in Table [2] 

This second step allowed us to select possible satellites 
and reject all other background or foreground objects charac- 
terized by a radial velocity difference (with respect to the pri- 
mary galaxy) more than V ou , = 1000 km s" 1 (McKay et al. 2002; 
B05). We refer to these galaxies as candidate satellites through- 
out the paper. 

In designing the masks containing the slits, we did not con- 
strain the wavelength range. Slit positioned on the north side 
of each quadrant produced spectra in the red range, with wave- 
lengths longer than 5500 A. Slit positioned on the south side 
of each quadrant produced spectra in the blue range, with wave- 
lengths shorter than 6500 A. Our strategy is to derive the redshift 
of satellites either from the Ha, [N II] emission lines or from the 
Hj3, [O III] emission lines that always fall in the observed spec- 
tral range. We were typically able to identify about 150 redshifts 
for each primary galaxy. 
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Table 2. Emission and absorption lines used in the preliminary 
redshift determination. 



Line 


Wavelength (A) 


[Oil] 


3726.03 


[Oil] 


3728.82 


CallK 


3933.66 


CallH 


3968.47 


H5 


4101.73 


Hy 


4340.36 


Hj3 


4861.31 


[Oiii] 


4958.92 


[Oin] 


5006.84 


Mg 


5183.6 


Ha 


6562.80 


[Nil] 


6583.41 


[Sll] 


6716.47 


[Sn] 


6730.85 



After the satellite candidates of a primary galaxy were iden- 
tified, we proceeded with the accurate measurement of their ra- 
dial velocity. We again, used all available emission lines (typi- 
cally Ha, [Nil] or Hj3, [Oil]) and used the sky emission lines 
to refine the wavelength calibration. We first spatially averaged 
the spectrum of each satellite and derived its one-dimensional 
spectrum. We then measured the wavelengths of a number of 
night-sky and target emission lines. We later, computed the dif- 
ference between the expected night-sky wavelength values and 
the measured ones, which were fit with a one-degree polyno- 
mial. This was used to correct the wavelength of the satellite 
emission lines. To have independent measurements, when possi- 
ble we measured the radial velocity separately for each exposure 
and for each emission line and then averaged the measurements. 
The radial velocities were measured with a typical uncertainty 
of 10 km s _1 . All radial velocities were reported with respect to 
the heliocentric rest-frame. 

The last step consisted of identifying the positions of all the 
200 slits on both the pre-imaging frame and the SDSS images 
and, in particular, the positions of the satellites. 

In view of future observations with the same instrument, we 
note, that we can easily increase the fraction of observed extra- 
galactic sources, and therefore satellites, by using different MOS 
masks for the same primary instead of repeating the same mask 
several times as we did in this pilot project. In very few cases, 
we added the spectra before measuring the redshift and we found 
that two independent 45 min spectra per target would have pro- 
vided a velocity measurements with uncertainties not significa- 
tively larger than what we obtained with several exposures. 

3. Distribution of satellites 

At the end of the data reduction phase, we found a total of 77 
objects (considering both VIMOS and B05 data) that fall in a 
projected circle in the sky of about 1000 kpc of radius and have 
a velocity that is at most 1000 km s smaller or larger than 
that of the primary galaxy. In Table [3] we indicate the num- 
ber of satellite candidates found for each primary galaxy (for 
J003828+000720 the quality of the data was insufficient for 
the analysis). Not all objects are physical satellites since out- 
liers may be present (see Sec l3.2l >. To assess the completeness 
of our sample, we estimate the fraction of objects that we ob- 
served over the total number of objects present in the field. For 
each VIMOS field, we derived the r-band magnitude distribu- 
tion of all extended sources, regardless of their redshift, from 



the SDSS DR7 photometric catalog. We then identified in the 
latter the objects that we had observed and found that, in the 
magnitude range between 18 and 21 r-band, they are about half 
of those present in the SDSS photometric catalog (see Fig. O 
and this ratio is roughly constant with the apparent magnitude. 
Only for objects fainter than 21 the ratio start to increase. The 
discontinuity visible at magnitude -~ 16 — 17 is caused by our 
having for these objects an SDSS redshift and we selected only 
the few galaxies with a recession velocity similar to the one 
of the primary galaxy. We can assume that, by adopting differ- 
ent masks in our observations, we could spectroscopically ob- 
serve twice the number of objects we actually observed. Since 
we placed the slits on targets with no particular criteria other 
than that they are extended sources, the real number of candi- 
date satellites (galaxies with a radial velocity V r in the interval 
V primary — 1000 km s _1 <V t < V primary + 1000 km s -1 ) present 
in the field is about twice the number we counted. We have to 
take into account this factor of about two when estimating the 
number density of objects. 

Table 3. Number of candidate and confirmed satellites found for 
each galaxy. 



ID 


SDSS a 


SDSS+VIMOS b 


Satellites' 1 


J134215+015126 


4 


8 


7 


J14521 1+044053 


4 


7 


7 


J152621+035002 


4 


17 


16 


J153221-002549 


5 


19 


5 


J154040.5-000933 


7 


11 


11 


J 154904-004023 


2 


8 


8 


J221957-073958 


4 


7 


7 



number of candidate satellites known from SDSS (B05). 
k candidate satellites from SDSS plus VIMOS observations. 
c satellites after rejection of interlopers. 



3.1. Structural properties of the primaries from their rotation 
curves 

A crucial advantage of this study is that by measuring the 
RCs of the primaries, we can estimate their inner gravitational 
field. As a first step, from their velocity-position diagrams 
and inclinations i, we derive their corrected rotation curves. 
Secondly, we compute the average velocity of all curves at ra- 
dius R opt = 3.2 Rjy (where Rq is the exponential disk scale- 
length). We define this velocity (V(R opt )), and we use it to 
normalize every RC to a same representative amplitude and ra- 
dius. We found (V{R opt )) = (255 ± 16) kms" 1 and (R opt ) = 
(19.8 ± 2.0) kpc. Because of the small luminosity range of the 
primaries, both quantities lie in a relatively small range of values, 
in agreement with the Tully-Fisher relation. Following Persic, 
Salucci & Stel ( 1996; hereafter PSS), for each primary we plot 
V(R(R opt )/R pt) (V(R op t)) /V(R ov t) (see Figs. E] and Q3, a 
quantity which is proper when comparing rotation curves of ob- 
jects of (slightly) different luminosity. We find, in agreement 
with PSS and Salucci et al. ( 120071 hereafter S07), that the op- 
tical RCs of the primaries have similar shapes and define uni- 
versal curve Vunc{R/Ro-puL) that accurately represents each of 
them. This implies that these galaxies have similar disks em- 
bedded in similar dark matter haloes. Therefore their RCs and 
the recessional velocities of their satellites probe the dark matter 
halo around a spiral of luminosity (Mr) with characteristic ra- 
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Figure 2. Lower panel: SDSS DR7 number counts of galaxies 
averaged in a 14' x 14' field around each primary galaxy {blue 
dotted line), number counts of the galaxies that have been ob- 
served with VIMOS (black line, magnitudes from SDSS), and 
number counts of the physical satellites (shaded region). Upper 
panel: ratio of the number of galaxies present in the field to the 
one observed with VIMOS. 

dius (/? pt)- The virial mass, virial radius, and halo escape veloc- 
ity can roughly be estimated by extrapolating its representative 
RC out to large radii according to V(R) = V(R opt ){R/R opt y - 2 , 
i.e. by adopting the velocity profile suggested by observations in 
S07 and found in ACDM simulated haloes. 

We find M vir - 3.5 x 1O 12 M and R vir ~ 400 kpc, which we 
consider as indicative values of the virial masses and radii of the 
primaries. Furthermore, we can estimate the escape velocity as 
V e (r) = f™V(r) 2 /r dr, because it mainly depends on the gravi- 
tational potential at small radii, which can be accurately mapped 
by the (co-added) RC. This estimate is robust in spite of the un- 
certainty in the RC extrapolation: even if we truncated the halo 
at 1/2 R v i r , V e (r) would be smaller by only 20%. We also note 

1 /3 

that the value of V e (r) varies with halo mass only as Mj ir . 

We use these estimates just as templates of these quantities 
and V e (r) helps us to distinguish interlopers from true satellites. 
We note that the actual circular velocity curve leading to the 
gravitational potential for our galaxies will be derived from the 
velocity dispersion of the satellite population. No result of this 
paper depends on the details of the above RC extrapolation. 

3.2. Interlopers 

Not all the galaxies detected around the seven primaries are their 
satellites: the distribution of objects around a target galaxy is 
contaminated by the presence of a background population of 
interlopers/outliers, i.e. of objects physically unrelated to the 
galaxy targets but at a small projected distance from them. That 
we investigate haloes of approximately the same mass and virial 
radius, and that these quantities are roughly known, allows us 
to infer the membership status of nearby objects according to 
their location and relative velocity with respect to the primary. 
Moreover, we can also distinguish the background of "false 
satellites" in outliers, galaxies around the primaries whose large 
value of relative velocity V re i show that they are systems un- 
bound to the primaries, and true velocity interlopers, i.e. objects 
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Figure 3. Satellites velocities with respect to the primary (open 
circles) versus radius. Black points are outliers. The red curves 
are the caustics identifying satellites and outliers. Vertical red 
lines indicate the bin locations and sizes. 



with V re i and projected locations similar to those of true satel- 
lites. 

In general, the removal of interlopers is quite complicated 
and uncertain. However, according to Chen et al. (120061) and 
Azzaro et al. (2006) for haloes of 10 12 - 10 13 Mq, the removal 
of outliers can be routinely done and the contamination of ve- 
locity interlopers becomes important only at projected radius of 
R> 400 kpc - R vir . 

A modern and efficient way of identifying member satellites 
is the caustic method introduced by Diaferio (1999) (later im- 
proved in Serra et al. ( 120101 1201 11 1). The method uses relative 
velocities V re i and projected radii to delineate the caustic curves 
that represent the escape velocity of the system (for details, see 
Serra et al. (120101 1201 IV ). We applied it to our data identifying 
the interloper galaxies and removed them from our sample. In 
Fig. [3] we present the V re i versus R for all objects, along with 
the caustic defined by Diaferio ( 1999). We also show the loca- 
tions of the bins adopted (ten objects per bin) in computing the 
satellites mean velocity and velocity dispersion. 

Since the present work is limited by the relatively small num- 
ber of satellites, we investigate this issue by means of a further 
method to identify objects gravitationally unbound to our pri- 
maries. We also used the negative binding energy criterium: we 
resorted to the negative binding energy criterium that uses the 
V e (r) profiles by obtaining a very similar result in the selection. 

With the above considerations, we removed from our sam- 
ple, at any projected radius R, the unbound/unrelated objects to 
the primaries whose inclusion strongly affects the estimate of the 
velocity dispersion. These outliers are galaxies more distant than 
5 Mpc from the primaries, which happen to lie on a line of sight 
passing close to their centers (i.e. < 1000 kpc). In Fig. [4] we 
give the histogram of our objects after we performed the velocity 
clipping described above, i.e. we show the number of candidate 
satellites for each radial bin and the number of outliers. Because 
of small numbers, the cosmic variance of outliers cannot be es- 
timated in this study. 

The clipping does not remove objects that are physically 
unrelated to the primaries, but reside at a smaller distance (1 
Mpc < D < 5 Mpc). In this case, V re i, which is the resultant 
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of redshift and peculiar velocity, can take similar values in a 
true satellite or in a physically unrelated object. This popula- 
tion of velocity interlopers affects in particular the satellite's sur- 
face density (at outer radii); we investigate its bias of our mea- 
sured projected satellite surface density by following Azzaro et 
al. (2006 ) and assuming that it provides a constant background 
(Chen et al. [20061 Azzaro et al. [20061 Mamon et al. l20T0l >. 
We estimate the latter with the SDSS r-band luminosity func- 
tion (Blanton et al. 120031 1. i.e. a Schechter function 4>(M) = 

0. 4<&*log(10) io - 4 ( M *- M )( a + 1 )exp(-10°- 4 ( M *- M )), with the pa- 
rameters <I>* = 0.0074 Mpc~ 3 , M* = -21.31, and a = -0.89. 

In detail, we estimate the number of velocity interlopers 
N V i(R), i.e. of the objects that lie within a projected distance R 
from the center of a primary galaxy but at a redshift distance 
D e < Vf/72 Mpc. Assuming that V e is constant with radius, we 
get 

Nvi = Pil 5 ^(M)dM / S 2%V e {R' /12)R'dR' 
= 1 x (fl/400kpc) 2 , 

where Mr = —14.5 is the limiting magnitude in our observa- 
tions, and 72 km s~ ! Mpc is the Hubble constant. In agree- 
ment with Azzaro et al. (2006]!, we consider the radial depen- 
dence of V e (see Sec.[3j} and find that/V,,, ~ 1 X (R/l Mpc) 175 

1. e. in our case the number of interlopers inside R v i r ~ 400 kpc is 
at most one, while in the region between Z? VI > and 2R v i r (that cor- 
responds to the three outermost radial bins) the expected number 
of interlopers is a few. In terms of number surface density, we 
find that jit,.,- ~5x 10~ 7 kpc -2 , with an uncertainty of a factor 
of two. Thus, the number of objects with Mr > —14.5 after a 
V e clipping represents a fair estimate of the number of satellites 
around spirals of Mr ~ —22. 



The first can be written as 

log M(*)-log Hm(R) = -4.44 + 0.18 x (*/100) + 

-7.0 x 10- 2 (fl/100) 2 ' 1 ' 

where R is in kpc and [I in kpc -2 , the above holding inside ~ 
500 kpc. We can remove the velocity interlopers, although they 
contribute in a negligible way for R < 700 kpc thus, to a first 
approximation, we take Eq. (|2]i as the estimate of the surface 
density of satellites. 

The projected radial distribution of galactic satellites around 
high luminosity primaries has been measured with the SDSS 
(C hen et al. l2006l l and the 2dFGRS (Sales et al. 120051 Sales et 
al. (|2005) also distinguished the case of blue (spiral) and red (el- 
liptical) primaries. In both surveys, the limiting flux of the satel- 
lites is 2-3 magnitudes brighter than that of our study and halo 
masses around primaries are a factor 2-5 smaller than the mass of 
the halo around our primaries. In view of this, the smaller ampli- 
tude of the satellite surface density found in the quoted surveys 
is unsurprising. It may indicate a bias in the radial distribution of 
satellites that depends on their luminosity. We see that fi (R) is 
not a power law: there is a central region of constant value, ex- 
tending for ~ 100 kpc. A similar trend can be seen in the quoted 
surveys, especially for the subsamples of red primaries in 2dF- 
GRS, and high luminosity primaries in SDSS. 

Outside 500 kpc, fl(R) is subject to the increasingly limited 
quality of the statistics. However, the primary inner kinematics 
guarantees that the density profile represents objects related to 
the primary but not virialized. The distribution of satellites as- 
sociated with an 3 x 10 12 M© halo seems to be cut off at about 
700 kpc. The observed fall-off in number density (see Fig. [5j is 
real: we could have detected objects at R > 2R vlr , if they existed. 




200 400 600 800 
R [kpc] 



1000 



Figure 4. Number of candidate satellites (top panel) and outliers 
(bottom panel) as a function of radius. 



3.3. The radial distribution of satellites 

The radial distribution of satellites was analyzed with bin sizes 
of ~ 100 — 200 kpc (ten satellites per bin), starting with the to- 
tal number of satellites detected around all primaries. To get a 
representative number for an individual galaxy, the frequencies 
were later multiplied by the factor of 2/7, where 7 is the number 
of primaries in the present work and 2 is the completeness factor 
discussed above. In the full project, we aim to reduce the latter 
factor down to ~ 1. In Fig. [5] we plot the number surface den- 
sity jJ.(R) before (filled circles) and after (full line) the removal 
of interlopers. 




Figure 5. Surface density of candidate satellites (filled circles) as 
a function of radius. The full line is the surface density of satel- 
lites after accounting for position interlopers. The dotted line 
indicates the surface density radial slope of a NFW halo. The 
full and open triangles represent the surface density of satellites 
found by Sales et al. (2005) and Chen et al. (2006), respectively. 
The dashed lines show the estimate of the interlopers velocity 
for /I = 3 x 10~ 7 kpc -2 and n = 10~ 6 kpc~ 2 . 
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Figure 6. Number of satellites against position angle differences. 
The total distribution is given by the upper histogram, while that 
of approaching satellites is represented by the unshaded area. In 
each bin, the rest of the contribution is due to receding satellites 
(shaded area). 



The satellite surface density distribution cannot be repro- 
duced by a projected Navarro-Frenk- White (NFW) profile (jj. 
R~ 1 - 5 ) for 0.1 < R/R V ir < 1 (see Fig.[5]l. In other words, it is not 
compatible with the predicted density distribution of dark matter 
haloes in a ACDM scenario. In particular, the predicted ~ 200 
subhaloes around a 3 x 10 12 M© halo exceeds the number of ob- 
served satellites by one order of magnitude, a situation we also 
find for Local Group galaxies. However, astrophysical processes 
in the late stages of galaxy assembly might depopulate most sub- 
haloes (e.g., Maccio et al. 1201 Oi l and the distribution of visible 
satellites might have no resemblance with that of existing DM 
sub-unities. 

We might try and compare our radial distribution with those 
predicted by state-of-the-art simulations of galaxy formation 
with models based on cold dark matter (CDM). The results of 
simulations using N-body techniques as well as different semi- 
analytic models of galaxy formation can be found in the litera- 
ture (e.g., Moore 1999; Baugh 2006). However, these codes do 
not produce unique results. We can add as a general comment, 
that these models may have difficulties in explaining the rapid 
increase in the cumulative number of satellites inside 100 kpc 
(Maccio et al J20101 >. 

3.4. Angular distribution of the satellites 

From the /?-band photometry, we measured the position angle 
(PA) of the primaries and derived APA = PA sate m te — PAga^y 
for each satellite. In Fig. [6] we plot the number of satellites as a 
function of APA and in Fig. [7] APA as a function of the distance 
to the primaries. We compared the PA distribution to a uniform 
distribution by means of a Kolmogorov-Smirnov test. The prob- 
ability that the PA distribution and the uniform distribution have 
a common origin is ~ 98 % (see Fig. [8). Hence, the observed 
distribution is consistent with being isotropic at current levels of 
uncertainty. In any case, it is clear that the distribution of satel- 
lites around primary galaxy in Figs. [6] and [7] do not show any 
preferential alignment, i. e. they are distributed isotropically. 
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Figure 7. Distribution of satellites around primary galaxies. 
Radius and position angle with respect to the major axis are 
shown for approaching [open circles) and receding (crosses) 
satellites. 
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Figure 8. The cumulate vector of APA (broken curve) is com- 
pared to that of a uniform distribution (straight line). A 
Kolmogorov-Smirnov test shows that the probability that the two 
distributions are coming from and identical underlying one is 
-98%. 



3.5. Velocity dispersion profile and rotation velocity at large 
radii 



Table 4. Velocities and dispersion versus radius. 
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<v re i) b 


o c 


kpc 
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521 


101 


32 
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-72 
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799 


86 
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average radius of the bin 
k average velocity of the bin and its uncertainty 
c velocity dispersion o~ of the bin and its uncertainty 
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Figure 9. Mean velocity of satellites as a function of radius. 
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Figure 10. Velocity dispersion as a function of radius. 



dure would of course also smear out any velocity anisotropies 
if they had different structure galaxy by galaxy. In this frame- 
work, we derive the circular velocity V = Rd<t>/dR where <t> is 
the gravitational potential, from the velocity dispersion in two 
different ways. 

First, we apply the usual argument that velocity dispersion 
has three degrees of freedom while rotation has just two, and 
we set V(R) = {'i/2) l l 2 a(R). Rotation velocities computed this 
way are plotted in Fig. QT| together with the inner RCs. 

Second, we adopt the fully correct procedure: we solve the 
Jeans equation for the 3D distribution of satellite n s (R) that we 
get from an Abel-integral deconvolution of )J.(R). We note that 
in this pilot study the quality of the statistics is still low and this 
is done just to illustrate the feasibility of the full project (which 
involves several times more data). The deconvolution yields 



V(R)=k(R)a(R) 



(4) 



with k = {-dlog n s {R)/dlog R-2 dlog a /dlog R) 05 . For r < 
400 kpc, i.e. for velocity dispersions of the three innermost bins 
we find k — 1.0,1. 4, 1.8. For the outer three bins in which very 
likely the Jeans equation cannot be applied, we allow V to range 
between (3/2) 1 ' 2 a and a (an object radially infalling) and we 
do not estimate a value. Fig. [T2] shows the inner RCs together 
with velocities computed with this second procedure. 

The last two bins illustrate a steep, large increase in the ve- 
locity dispersion. This is likely to indicate that most, perhaps 
all objects, are not physically bound to the primary galaxy. The 
detected sharp decrease in the number of true satellites might 
indicate at 600 kpc the edge of the dark matter halo. With the 
full sample that is four times larger, we hope of course to have 
consequently a more accurate determination of the mass distri- 
bution. 



We estimate the velocity dispersion o(R) of the seven co- 
added systems by binning the satellites into groups of ten ob- 
jects. For each bin, we compute the average velocity (V re i) and 
the velocity dispersion a. We report them in Table [4] and plot 
them in Figs. l9l and [TUl The average difference (V re i) is compat- 
ible with zero at any radius (see Fig. [9]). In addition, along the 
line-of-sight the velocities are isotropic. This is confirmed by the 
finding that (V re i) ~ 0.2(7 at any radius (see Table|4]l. The flat ra- 
dial trend also implies that the satellites system has a negligible 
rotation. 

In view of the considerably large r.m.s. of cr(R), we chose 
the simplest linear function to fit data from 100 kpc to 500 kpc 

a(R) = [265 -0.35(7? -100)] km s _1 , (3) 

where R is in kpc. 

The inner RCs of the primaries indicate, following Salucci 
et al. (2007), for our co-added system a virial mass of ~ 3.5 x 
10 12 M and a virial radius of 400 kpc. Inside this radius, <7 
measures the kinetic energy of a unit mass test particle in grav- 
itational equilibrium inside the galaxy potential. At larger radii 
in the range 400 < R < 800, kpc we have two cases: 1) a still 
measures the kinetic energy, but of objects out of equilibrium or 
2) cr is totally contaminated by velocity interlopers. 

The velocity dispersion measures the circular velocities well 
beyond the extent of the rotation curves. We assume that the 
satellites have isotropic motions. This assumption is supported 
(albeit not proven) by the isotropic angular distribution and ve- 
locity dispersion along the line of sight. The co-addition proce- 



3.6. Mass modeling 

The radial profiles of circular velocities shown in Figs. [TT] and 
Q2]were fit with the universal rotation curve (URC) with (i?rj) = 
6.2 kpc and halo mass M vir = 3.5 x 10 12 M that closely fits 
the inner RCs. The URC has no free parameters to fit the present 
outer data, therefore its success is remarkable. The URC depends 
only on the disk scale-length Tfo and the halo mass M v ; r , which 
also controls the disk mass (Shankar et al. 2006). The URC mass 
model fits well the data points in both cases for M,,„ = 3.5 x 
10 12 M R . 



4. Discussion and conclusions 

In this pilot project, we have searched for satellites around seven 
spiral galaxies of comparable luminosity (Mr ~ —22). We have 
found 61 objects that, based on their velocity relative to the pri- 
mary galaxies, are their satellites. By combining results, we have 
investigated several properties of these satellites, which allow us 
to draw a number of conclusions. 

We find that our target dark matter haloes have about 15 
satellites (down to Mr = —16) inside their virial radius. Their 
measured inner kinematics show that, as intended, they have 
similar halo masses, of about 3.5 x 10 12 M and the embedded 
luminous spirals have luminosities Mr ~ —22. Therefore this 
result is comparable to what we observe for M3 1 and the Milky 
Way (MW) in the Local Group. As in those cases, this small 
number of satellites represents a problem for "naive" ACDM 
theories, which predict an order of magnitude larger number of 
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Figure 11. RCs of primaries and circular velocities obtained 
from the velocity dispersions of the satellites (black squares) 
(V(R) = (3/2) 1 / 2 a(/?), see text for details). The URC is shown 
(as a black line) out to the virial radius. 
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Figure 12. RCs of primaries and circular velocities obtained 
from velocity dispersions and jeans Equation (black squares, see 
text for details). The URC is shown (as a black line) out to the 
virial radius. 

objects orbiting around dark haloes of this size. To reconcile ob- 
servations and theory, one has to assume that most of the galaxy 
subhaloes have become unobservable. This possibility, related 
also to the apparent lack of subhaloes around the Milky Way, 
is being investigated (e.g. Kravtsov et al. 12004; Libeskind et al. 
2009). Could a long list of physical processes occurring in sub- 
haloes during their history jeopardize (while perhaps also intro- 
ducing a radial bias) their ability to host luminous objects? In 
this respect, we note that our primaries have a mass 3-4 times 
higher that of the MW and the combined sample has already 10 2 
objects and in the future will have ten times more. We will be 
able to investigate the "missing subhaloes problem" also in sub- 
haloes with a mass 1 order of magnitude larger than that of the 
Large Magellanic Cloud (LMC), and therefore in subhaloes that 
should host normal and not dwarf galaxies. It will be interesting 
to see if, at these high masses, the problem remains, and can be 
tackled by the explanations put forward in the case of MW. 

By considering a composite system including all detected 
satellites, we are able to study their radial distribution with re- 



spect to the central galaxy. The surface density has a shallower 
profile than found in previous surveys, based on fewer satellites 
per primary. We propose that this is due to a bias in the radial 
distribution of satellites that depends on their luminosity. 

We also find that the angular distribution of satellites around 
primary spirals is isotropic. Furthermore, the relative velocity 
along the line of sight is isotropically distributed around zero. 
Based on this observation, we can use the velocity dispersion 
profile to estimate rotation velocities far beyond the extent of 
the inner rotation curves. 

An additional result of this study is that the RCs of the pri- 
maries appear very similar. Because we targeted objects with the 
same absolute luminosity, our finding confirms that RCs mainly 
depend on an object's mass. It follows that the mass distribution 
inside the optical region is the same for all our objects. It also 
justifies combining the satellite relative motions to derive the ve- 
locity dispersion and probe the mass distribution of a "stacked" 
halo. The outer RC obtained from the velocity dispersion results 
is smooth and continuous with respect to the inner measured RCs 
and, over two orders of magnitude in distance, is very well re- 
produced by the S07 universal mass model that holds for spirals. 
In this paper, the model is tested for the first time with kine- 
matical measurements out to the virial radius of the galaxy. The 
technique we have applied is therefore very promising. 

In the future, we plan to test our findings by enlarging our 
sample of primary galaxies. We estimate that by co-adding the 
satellite distribution and kinematics of ~ 20 systems of a same 
high mass we can infer the mass distribution of the underly- 
ing DM halo. The full project aims to replicate this procedure 
for haloes of intermediate and low mass (8 x 10 11 M and 
2 x 10 M respectively), to derive the mass dependence of the 
DM distribution. 
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